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IMP I OPTICAL ASPECT SYSTEM 


by 

E. John Pyle 

Goddard Space Flight Center 


INTRODUCTION 

The problem of attitude or aspect determination is the problem of defining the angular relation- 
ship between two coordinate systems, one fixed in a reference space and the other fixed in the space- 
craft. There are six degrees of freedom associated with a general coordinate system, with respect to a 
reference coordinate system. Three of these concern the translational motion of the center of mass of 
the spacecraft, which is a problem in orbit and trajectory calculations. The remaining three degrees of 
freedom deal with the rotational motion of the spacecraft about its center of mass. The measurement 
of this rotational motion is called “aspect determination.” 

A body spinning in free space has most of the properties of a gyroscope; if a spin is imparted to a 
spacecraft at the time of injection into orbit, it performs like a gyroscope while traveling in free space. 
An aspect system included in the spacecraft assembly gives the orientation of the spacecraft’s spin axis. 
Since there are three degrees of rotational freedom of a rigid body associated with aspect determina- 
tion, the instrumentation in the aspect system must be capable of measuring at least three parameters. 

This paper discusses the IMP I Aspect System; Figure 1 is an overall block diagram. The param- 
eters measured are the elevation angles— with respect to the spin axis— of the sun, )3, and the earth, 5; 
the order of observation; and the fraction of one spin period between the observation of each. Addi- 
tional information such as the position of the spacecraft with respect to the two observed references is 
required ; this is available from earth-based observations. 

The order of observation and the fraction of a spin period between observations can be measured 
on board the spacecraft. These may be derived from a knowledge of the exact times that a specified 
reference plane crosses the solar and terrestrial disks. These times must be related to some known time 
standard; these relationships must be remembered, encoded, and presented to the transmitter. Two 
sensors will be used to determine the angle between the spin axis of the vehicle and each of the two 
references. Since the sun, for practical purposes, is a point source to the sun sensor, the vector from 
the spacecraft to the sun is an axis of symmetry. This means that the locus of the spin axis is a cone 
whose axis is the sun vector and whose half-angle )3 is the angle measured by the sun sensor, that is, the 
angle between the spin axis and the sun vector. In order to calculate the position of the spacecraft’s 
spin axis, a measurement must be made with regard to another reference— the earth. The aspect system 
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Figure 1 -System block diagram. 





therefore will contain an earth sensor, sensitive to visible light, with a narrow field of view mounted so 
as to scan 90 degrees from the spin axis. Information from this sensor will lead to the calculation of 
the angle between the spacecraft spin axis and the vector from the center of the earth through the cen- 
ter of gravity of the spacecraft. This vector, also, is an axis of symmetry, and the locus of the spin axis 
is a cone whose axis is this earth-spacecraft vector and whose half-angle 6 is the angle calculated from 
earth-sensor information. Hence the spin axis must be at one of the points of intersection of the two 
cones (see Figure 2). An unambiguous solution is obtained when the two other parameters measured 
on board the spacecraft are taken into account. 

The Spin Synchronous Clock (SSC) is a digital system that generates a time-based clock signal 
synchronized to the spin rate of the spin-stabilized spacecraft. The basic function of the device is to 
generate 2" (where n = 7) pulses per spin period. These pulses divide the spin period into 2” - 1 equal 
time intervals plus one remaining interval, which may be slightly greater or smaller than the other in- 
tervals. To accomplish this function the SSC receives a high-frequency signal f Q = 25.6 kHz from the 
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Figure 2— Relationship between the spacecraft spin axis and the intersection of the two 

locus cones. 



spacecraft encoder. The command pulse defining the rotational period is generated in the optical- 
aspect electronics and fed to the SSC. 

ASPECT SENSORS 

The energy emanating from a celestial body can be detected by sensitive measuring instruments. 
The optics, fields of view, and other characteristic parameters of these sensors are predicated on the 
energy source selected and the particular sensor application. 

Attitude can be derived from the spin-axis orientation relative to any number of celestial objects. 
A minimum of two objects (sun and earth) is necessary to completely specify the spin-axis attitude of 
the spacecraft. The sun can be treated as a point source, whereas the earth’s angular subtense varies as 
a function of relative sun position and closeness of the detecting instrument. 

Digital Solar Sensor 

The digital solar sensor (Figure 3) measures the angle of incident sunlight with respect to the sen- 
sor Z-axis and expresses this angle as a digital number. The incident sunlight, passing through a slit on 
the top of a quartz block, is screened by a gray-coded pattern on the bottom of the block to illuminate 
or not illuminate each of the photocell detectors. The angle of incidence determines which combina- 
tion of photocells is illuminated. The solar sensor also includes a command slit, which is mounted 
perpendicular to the gray-coded reticle. If the sensor is rotated about a vertical axis along the com- 
mand slit, the field of view of the two slits will sweep over a solid angle. When the plane containing 
the command slit passes across the solar disk, one or more of the photocells will be illuminated. The 
time that this illumination occurs provides a measurement of the azimuth angle of the sun. The par- 
ticular combination of photocells that is illuminated provides a digital measurement of the elevation 
angle of the sun in sensor coordinates. 

The digital solar sensor has features that make it more desirable than the analog types of solar 
sensors: (1) It is not subject to errors introduced by earth shine, (2) there are no components that can 
drift, (3) no in-flight calibration is required, (4) the weight is small, and (5) power is needed only to 
drive the output load. 

Earth-Horizon Detector 

The horizon detector consists of a simple telescope and lens. The detector element is a photo- 
diode, which is placed at the focal point of the lens. The field of view of the horizon detector is a pen- 
cil beam approximately 1 degree in diameter. If the detector is mounted at an angle, y, from the spin 
axis of a rotating spacecraft, the field of view of the detector will traverse the surface of a cone whose 
half angle is equal to y. When the field of view scans the discontinuity caused by the sunlit earth 
against the dark background of interplanetary space, this detector generates an output signal. This sig- 
nal, after amplification, is differentiated to form pulses at its leading and trailing edges. The informa- 
tion contained in the relative position of these pulses to the command slit’s crossing of the sun provides 
a measure of the inclination of the spacecraft zenith vector to the spin axis. 
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Figure 3— Schematic representation of a digital solar sensor. 


ASPECT-SYSTEM ELECTRONICS 

Solar-Sensor Electronics 

The solar-sensor output consists of 10 channels of sun information. Nine of these contain angular 
information; the 1 0th contains time-occurrence or azimuth information. 

Amplification of the output signals of the nine angular information channels produces a positive 
voltage when the channel is excited, i.e., when sunlight is incident on the detector (see Figure 4). 

These amplified signals are then placed in storage and, at the proper time, shifted through an inverter 
to the spacecraft encoder. 

The amplified azimuth information, referred to as the “solar command pulse,” occurs when the 
sun is in the field of view of the sensor command slit. The width of this pulse is a function of the spin 
rate and the command-slit field of view. This is undersirable; to eliminate these effects on pulse width, 
a centered sun pulse is generated at the midpoint of the solar command pulse (sun-centering module). 








The logic that generates the centered sun pulse also produces the transfer (load) pulse, which stores the 
sensor’s angular information (aspect module C). 

The centered sun pulse is used to measure sun time and spin period. Sun time is the interval be- 
tween a known spacecraft encoder function (C24) and the following centered sun pulse. Spin period 
is defined as the time interval between two successive centered sun pulses. 

Earth- Detector Electronics 

The earth detector is a light-sensing device that detects the gradient between the sunlit earth and 
outer space. The sensing element is a photodiode having a response that peaks at 1 micron. As the 
spacecraft rotates, the detector’s field of view sweeps over a cone. At certain points in the orbit, the 
scan of this detector successively crosses the horizons and/or terminator of the sunlit earth. When this 
occurs, the detector produces an electrical signal. This signal is fed into a high-input-impedance ampli- 
fier, i.e., earth-amplifier module (see Figure 5). Subsequent electronic processing shapes this signal 
into positive pulses occurring at each intercept of the earth’s horizon and/or terminator. The same 
processing electronics eliminates any solar information that might be caused by the earth detector’s 
viewing the sun. 

The information contained in the relative position of these earth pulses to the centered sun pulse 
is used to compute the spin-axis earth angle, 5. The relative position of the earth and sun pulses are 
contained in two measureable quantities: earth time and earth width. Earth time is defined as the in- 
terval between the centered sun pulse and the earth pulse that defines the start of an earth scan. The 
second interval is the time between the two earth pulses, corresponding to the apparent width of the 
sunlit earth. 



Figure 5— Earth-amplifier module 16-175. 
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Measurement of Aspect Parameters 

The quantities measured by the aspect system are those already mentioned: 

(1) Spin period. 

(2) Sun time. 

(3) Earth time. 

(4) Earth width. 

(5) Spin-axis sun angle. 

As previously mentioned, the three parameters necessary for attitude determination are (1) j3, the 
spin-axis sun angle, (2) 8, the spin-axis earth angle, and (3) the fraction of a spin period between ob- 
serving each. 

Spin-axis sun angle (5 is measured directly by the solar sensor; its digital representation is placed in 
storage. 

Spin-axis earth angle 5 must be computed from the relative position of the sun and earth (see 
Appendix D). This information is contained in the five quantities measured by the aspect system. 

The fraction of a spin period between the observations of the two reference sources (sun and 
earth) can be determined from knowledge of the spin period, earth time, and earth width. 

Encoder Interface 

The spacecraft encoder’s interface with the aspect system places certain restrictions on the meas- 
uring, processing, and storing of aspect data. 

The encoder samples the aspect data once every C24 period, T C24 - 81.92 s. The encoder cir- 
cuitry supplies the storage accumulators for the interval measurements, i.e., sun time, spin period, 
earth time, and earth width. Since the nominal spin period for the satellite is 1 2 s, the gates to these 
accumulators must be controlled so that repeated data are not fed into them; i.e., only one measure- 
ment of a particular interval time is made in 81.92 s. 

To satisfy this timing restriction, the leading edge of C24 is differentiated and used as a synchro- 
nizing signal by the aspect system. The occurrence of this signal is used to initiate a series of events 
that provide only a single measurement of each of the above-mentioned time intervals in the 81.92-s 
encoder-controlled period. 

The gating of the sun-time accumulator is controlled by a flipflop that is reset by the centered 
sun pulse and set by the aspect synchronous signal. Hence, for the time interval between the aspect 
synchronous signal and the first centered sun pulse, this flipflop enables the accumulator and allows a 
1.6-kHz signal to be counted. Repeated occurrences of the centered sun pulse only serve to keep the 
flipflop in the inhibit state. 
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A two-stage counter, reset to zero by the aspect synchronous signal (aspect sync), controls the 
spin-period accumulator by generating an enable signal between the first and second centered sun 
pulses following reset. To prevent further generation of this enable signal during the C24 period, the 
input to this counter is inhibited following the second centered sun pulse. When the accumulator is 
enabled by its gate, a 1.6-kHz signal is used for a counting source. The trailing edge of this enable sig- 
nal is also used to control the transfer pulse which sets the angular information from the solar sensor 
into the storage register. 

The earth-time and earth-width accumulators are controlled by another two-stage counter. The 
first stage, which controls the earth-time accumulator, is reset by either the aspect sync or the earth’s 
leading-edge signal. The second stage is reset by either the aspect sync or the earth’s trailing-edge sig- 
nal and controls the earth-width accumulator. 

The action of this counter is initiated by a selected centered sun pulse. This pulse is generated 
once per C24 period by the first centered sun pulse following both the aspect sync signal and the earth 
amplifier signal. Hence, the data from the second earth scan in the C24 period is used to control the 
1.6-kHz counting signal for the earth-time and earth-width accumulators. 

The spin-axis sun-angle information is shifted out of its storage register once per C24 period by a 
series of shift pulses controlled by the encoder. The timing of these shift pulses is such that they occur 
well after the aspect system’s generation of the transfer (load) pulse. 

The above discussion was based on the encoder operation in its low bit rate mode (400 bps). The 
effect of the high bit rate mode (1600 bps) on the aspect data is simply to read out the same data four 
times in the C24 period. The added readout supplies no new data, since the C24 period generating the 
aspect sync signal does not change. 

Figure 6 shows in detail the timing for aspect parameters. 

SPIN SYNCHRONOUS CLOCK ELECTRONICS 

Logic Functions 

The basic function of the Spin Synchronous Clock (SSC) is to generate 2" (where n = 7) pulses 
per spin period. These pulses divide the rotational period into 2" - 1 equal time intervals plus one re- 
maining interval which may be slightly greater or smaller than the other intervals. To accomplish this 
function, the SSC receives a high-frequency signal f Q of 25.6 kHz from the spacecraft encoder and uses 
the centered sun pulse generated by the aspect system to define the rotational period of the spacecraft. 

The SSC includes the following components: a pulse generator; counters C Q , Cj, C 4 , C d , and C t ; 
a storage register; a comparator; and pseudo-command circuits. 

Generation of the basic control functions for the SSC is provided by two counters, C d and C 4 , 
and associated decoding gates. These basic control functions are as follows: 

(1) Cj inhibit signal— an inhibit voltage that prevents counter Cj from changing its contents dur- 
ing transfer. 
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Figure 6— Aspect parameters; timing diagram. 


(2) C 1 transfer pulse— a pulse that transfers the contents of counter Cj into a storage register. 

(3) Cj reset pulse— a pulse that resets counter C x to zero following the transfer pulse. 

The leading edge of the centered sun pulse (command pulse) resets both counters to zero, indicat- 
ing the start of a rotation period. Counter C d performs the frequency division on the f 0 signal to gen- 
erate the /j signal (/j = fj 2” = 25.6/128 = 0.2 kHz = 200 Hz). The block diagram, Figure 7, shows 
the logic for the three basic control functions and the fact that the input to C 4 is a gated function of 
one of the intermediate stages of C d . The gating insures that the basic control functions occur only 
once each spin period. Use of an intermediate stage of C d as the input to C 4 merely provides reset and 
transfer pulses of convenient duration. The only restriction is that these control signals be completed 
before the first count of enters the interval counter C r 

Upon the occurrence of the transfer pulse, the storage register accepts and retains the information 
from the interval counter C r Another counter, C 0 , along with the storage register, provides the inputs 
to the comparator. When the contents of the storage register and counter C Q are identical, the com- 
parator generates a pulse. In a given spin period, the pulses generated by the comparator become a 
clock signal synchronized to the spin period. These pulses, along with the command reset pulse (com- 
mand pulse), are used to reset counter C 0 . Hence, they must be such that their action does not cause 
the loss of an / 0 count. This condition determines the pulse widths of the comparator-generated 
pulses and the command reset pulse; i.e., the command reset pulse is the leading edge of the centered 
sun pulse. 
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The pulses generated by the comparator constitute the input signal to output counter C r Hence, 
decoding counter C t provides information as to the instantaneous pointing direction of a point on the 
spacecraft. 

During the rotation period, counter C t should reach a count of 2” (i.e., 128) and then be reset to 
zero by the command reset pulse. 

If a command reset pulse does not occur after 2" (i.e., 128) counts have been accumulated, the 
pseudo-command generator is activated, and the counter C t overflows and starts accumulating counts 
again. When C t next reaches 2 n (i.e., 128) counts and a command reset pulse has still not occurred, 
the pseudo-command generator produces a pseudo-command signal. This action continues, with a 
pseudo-command signal being generated every time C t reaches 2" (i.e., 128) counts. The occurrence 
of the command reset pulse causes C t to be reset to zero and the pseudo-command generator to be 
deactivated. Hence, the output counter C t will wait one spin period following the loss of the com- 
mand reset pulse before generating a pseudo-command signal for every 2" pulses from the comparator. 

System Performance 

The SSC divides the spin period into 2” sectors. To accomplish this, it uses the information in 
the storage register as the basic interval to be sectored. An investigation of the accuracy of this infor- 
mation will lead to an evaluation of system performance. 

Consider the representation of the spin period by interval counter C r This information may have 
a total error of minus one count of f v This follows directly from the division process performed in 
the generation of the signal. Since any division process can generate both a quotient and a re- 
mainder, the Cj representation can differ from the true spin period by the amount of this remainder. 
This means that if the spin period is not an integer multiple of T l = 1 jf v the contents of Cj will be 
less than the spin period. This difference can be as small as 1 or as great as 2" - 1 counts of f Q . 

Another consideration not previously mentioned is the interaction between the Cj inhibit signal 
and the f x signal when spin period is measured. This effect may cause an extra count to be placed in 
counter Cj if the inhibit signal occurs at a time when the f 1 signal is a “1”. This occurs when the con- 
tent of counter C d is greater than one-half its capacity,- i.e., when the last stage of C d is a “1”. 

The combined effect of these two considerations is to incorporate a round-off error into interval 
counter C r If the remainder of the division process is less than 1/2, the interval counter contains only 
the number of counts of f l that occurred during the spin period. If the remainder is equal to or greater 
than 1/2 of f v the interval counter contains a count greater by 1 than the number of counts of /j 
sampled during the spin period. Hence, the error in the interval counter has been translated from a 
count of -1 to a count of ±1/2 of /,. The total error in Cj as a function of the remainder in C d is 
shown in Figure 8. Thus when the remainder is less than 1/2 of f v the output counter attains a 2 n 
count before the next SSC command pulse. The maximum value of this time lapse is equal to TJ2. 

Another factor affecting system performance is the variation in two successive spin-period inter- 
vals. This variation may result either from a jitter effect of the onboard instrument’s generating the 
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Figure 8— Total error in C., as a function of the remainder in C d . 

command pulse or directly from the motion of the spacecraft (precession). Since the operation of the 
SSC is based on the assumption that two successive spin periods are equal, this spin-period fluctuation 
has a one-to-one effect on system performance. However, this effect simply reduces the accuracy of 
the location of the sector pulses in a given spin interval. 

A companion factor to the above short-term fluctuation is the long-term loss of the command 
pulse. This effect is taken into account in the operation of output counter C t . Counter C { operates in 
such a manner that, if the command pulse does not occur after 2" pulses, the counter will continue to 
count until either 2" +1 counts have been registered or the command pulse occurs. When 2” +1 counts 
have been accumulated, the associated circuitry generates a pseudo-command signal. This signal resets 
continuous counter C 0 but does not update interval counter C^. Therefore, the system waits one rota- 
tion period following the loss of the command pulse before generating a pseudo-command signal at the 
last known spin-period interval stored in the storage register. Hence, for a long-term loss of the com- 
mand pulse the system continues to operate, but with the last known spin-period measurement. 

Sector Determination: Location and Width 

In order to define the location in a given spin period of a specific sector region, a reference sys- 
tem must first be considered. Let the reference system be based on the 360 degrees of rotation in a 
spin period. (See Figure 9.) Then, each pulse generated by the comparator and defining the end of a 
sector region is separated from the previous one by 360/2" degrees of 'rotation. Hence, in this reference 
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Figure 9— Sector determination. 

system, the 2" sector pulse occurs coincident with the command pulse, and the width of each sector is 
360/2" degrees. 

The actual SSC generates sector pulses based on the spin-period representation present in the 
storage register. Hence, these sector pulses match the ideal reference system if, and only if, the re- 
mainder of the division process is equal to zero. In all other cases (where the remainder is greater than 
zero), the sector pulse generated by the comparator circuitry in the actual system is shifted from the 
reference system location by 


( T 0 Rk 
\ 2 " ’ 

) (2 n -R)T 0 k 
f+ 

V 2" 


for R < 2" _1 
for R > 2” -1 


where 

2" = number of sector pulses in a spin, 
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T 0 = 1 // 0 = period of the f Q signal, 

R = remainder (1 < R < 2" - 1), 

and 

k = number of the sector (1, 2, 3, ... 2” - 1, 2"). 

Now examine the case for R < 2" _1 , the case when the remainder is less than 1/2 of f v The stor- 
age register contains a time representation of the spin period which is less than the actual value. Hence, 
the interval being divided into 2" sectors is equal to the true spin period minus the time represented by 
the remainder (see storage register in Figure 7). Let 


Then, 


and 


N = number of counts in interval counter Cj . 


(2 n N + R)T q = true spin period , 

2” T q N = time representation of spin period by storage register , 


T a - 1 //„ 


Therefore, each sector pulse is generated slightly before the location of its corresponding ideal refer- 
ence sector pulse. This effect is accumulative and results in the 2”th sector pulse being generated T Q R 
seconds before the 2 n ideal reference sector pulse. Since each sector pulse resets continuous counter 
C 0 , this counter will have obtained a value of R counts when the command pulse occurs. 


Hence, if the first sector is defined as the time interval between the command pulse and the next 
pulse generated by the comparator, this results in an effective “dead time” between the end of the 
2”th sector and the start of the first sector in the next spin. That is, 


dead time width = RT n 


sector width = NT, 


o 


R < 2” _1 . 


Therefore, the time between the 2”th sector pulse and the first sector pulse in the next spin will in- 
clude the dead time and have a value of (TV + R)T Q . 


When R > 2”' 1 the information in the storage register represents the true spin period plus the 
effect of the inhibit signal adding a count to the interval counter C r This increases the contents of 
counter C x by 1 ; so N is replaced by N+ 1. 

We have 

(2 n N + R)T 0 = true spin period 

and 


2 n (N + 1 ) T 0 = time representation of the spin period by storage register . 


Hence, it is the quantity 2” (TV + 1)T 0 that is divided by 2”. This results in each sector pulse being gen- 
erated by the comparator slightly after its corresponding ideal reference sector pulse. Again this effect 
is accumulative and the 2”th sector pulse should occur a (2” -R)T Q time interval after its ideal refer- 
ence sector pulse. This time is also after the command pulse which resets both the continuous and 
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output counters. Hence the system will reset before the generation of the 2"th sector pulse. Hence, 
the comparator will generate 2” - 1 pulses in a given spin period. That is, 

dead time width = (2” - R)T n ) 

}r> r-K 

sector width = (TV + 1)7^ I 

Again, the dead time width is included between the 2" - 1 sector pulse and the first sector pulse in the 
next spin. This internal, referred to as the width of sector 1 , will have a value of [2(TV + 1 ) - (2” - T?))r o . 

Hence, the location of each sector pulse and the width of each sector can be determined from 
knowledge of the spin period and remainder values. 

Operating Range: / 0 , 2”, and Spin Ratew z 

The operating range of the SSC is determined by the relationship between / 0 in Hz, 2”, and the 
spin rate of the spacecraft co z in rpm. As previously mentioned, the division process can generate both 
a quotient and a remainder. The effect of the remainder is to shift the sector pulse from its reference 
location. The maximum shift occurs in the 2” sector and has a value of RT Q . Hence, if R is greater 
than the number of counts that represents the spin period in interval counter C 1? the shift caused by R 
may be greater than a sector width. Therefore, one constraint on the relationship between / 0 , 2", and 
spin rate is that the number of pulses in interval counter Cj be equal to or greater than remainder R. 
We have 


and 


JU.-z'-i 


iV. > 2” , 


where TV jc is the number of pulses fed into interval counter C r 

Since the number of pulses TV jc fed into the interval counter Cj is a result of frequency division 
by 2” , the minimum number of pulses TV jd fed into frequency-divider counter C d must satisfy 

TV > 2”TV 

■' v id L iv io ’ 

where TV, is the number of pulses out of counter C d . Therefore, 


^id = 


spin period _ 60/ 0 
r co 

o z 


^id - ^ic . 

N id >2 n N ic >2 2n , 

60 f 0 


CO 


> 2 ln , 
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and 


fo> 


2 2 «co z 
60 * 


A further constraint on the operating range of the SSC is that the number of pulses into the inter- 
val counter, N ic> must not exceed its capacity. Thus, 

2 k = capacity of interval counter C l (k= 17) , 

JV ic <2*, 


and 

and, since 
and 

we have 
and 


N io <2 k , 


N iH 

N- <— 
10 2" 


^id = ' 


60/ n 


CO 


60/ 0 t 

2 " -<2 k 

“z 

2 n+k u, 


fo<- 


60 


Consideration of both constraints yields 


2 2w co z 

60 


</ 0 < 


2 n+k o) z 
60 * 


System Output Signals 

The comparator generates the basic clock signal synchronized to spin rate. As previously men- 
tioned, this signal is subject to the effect of the remainder value, which causes noncoincidence of the 
2”th (i.e., 128th) pulse and the command reset pulse. Therefore, logic is provided to delete the 128th 
pulse from this signal and replace it with the command reset pulse. This results in the basic clock sig- 
nal’s having coincidence of the last pulse and the command reset pulse, independent of the remainder 
value. It is this coincident 128-pulse signal that is sent to the experimenter. (See Figures lOAand 10B.) 

The output of each stage of output counter C t is also available to all spacecraft systems. Separate 
buffering is supplied for each actual interface, and a special interface exists between the SSC and the 
spacecraft encoder. 
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Figure 10A— SSC output signals, basic logic. 
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Appendix A 


THEORETICAL CONSIDERATIONS OF ASPECT DETERMINATION 


Energy and Momentum Considerations 

In the solution to the problem of force-free motion of a rigid body, a point on the body can be 
defined in terms of the body-fixed axes, x, y, and z. For the purpose of defining the orientation of 
the body fixed axes relative to nonrotating reference axes, X, Y, and Z, Euler introduced three inde- 
pendent angles, 0, 0, and 8. These Eulerian angles are shown in Figure A-l and are defined by three 
successive rotations performed in a specific sequence. The sequence is started by (1) rotation of the 
initial set of axes X, Y, and Z, counterclockwise about the Z-axis through an angle 0, producing the 
intermediate set of axes, £', r?', and ?' (Figure A- la). (2) Next, this intermediate set of axes is rotated 
counterclockwise through an angle 8 about the £'-axis, producing a second intermediate set of axes, t 
rj, and J (Figure A-lb). The §'-axis is also called the line of nodes. (3) Finally, axes £, 17, and f are ro- 
tated counterclockwise about the f-axis through an angle 0, forming axes x, y, and z (Figure A-lc). 
Hence, the Eulerian angles, 0, 0, and 8, completely describe the orientation of the x,y, z coordinate 
system with respect to the X, Y, Z coordinate system. 

A mathematical method of representing the rotations described above and based on the matrix 
representation of an orthogonal transformation is given in Appendix B. 

Now let the X, Y, Z coordinate system define an inertial space. Assume a rotating spacecraft situ- 
ated in this inertial space where the momentum vector, L, of the spacecraft is orientated along the 
Z-axis of the defining coordinate system. Let the x,y,z coordinate system be aligned with the three 
principle moments of inertia of the spacecraft, I v I 2 , and / 3 . Hence, the x,y,z system is fixed in the 
rotating spacecraft. The instantaneous values of momentum about axes x, y, and z are, respectively, 
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and 


(A-l) 

(A-2) 

(A-3) 


Px =I : i°>x> 

Py =1 2°>y . 

Pz =/ 3 CJ Z > 

where a x , aj y , and co z are the respective instantaneous values of angular velocity about axes x, y, and 
z. The respective instantaneous values of the momentum p x , p y ,p z and the angular velocity oj x ,co y , 
cj z about the x, y, and z axes can be expressed in terms of the momentum vector (see Figure A-2 and 
Appendix C): 


and 


and 


P x -L sin 6 sin 0 , (A-4) 

P y ~ L sin 8 cos 0 , (A-5) 

p z = L cos 8 , (A-6) 

oj x = 4> sin sin 8 + 6 cos \p , (A-7) 

co y = 0 cos \J/ sin 8 - 8 sin \p , (A-8) 

co z = 0 cos 8 + 0 . (A-9) 

Substitution of Equations A-4 through A-9 in Equations A-l through A-3 yields 

p x = L sin 0 sin 8 = / x (0 sin 0 sin 8 + d cos 0) , (A-10) 

P y = L cos 0 sin 8 = / 2 (0 cos 0 sin 8 - & sin 0) , (A-l 1 ) 

p z = L cos 8 = / 3 (0 cos 8 + 0) . (A-l 2) 

Multiplying Equation A-10 by sin 0 and Equation A-l 1 by cos 0 and adding the results gives 

L sin 8 = (/j - / 2 )0 cos 0 sin 0 + 0 sin 8(1 x sin 2 0 + / 2 cos 2 0) . (A-l 3) 




Figure A-2— Components of angular momentum and velocities. 
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Multiplying Equation A- 10 by cos 0 and Equation A-l 1 by sin 0 and subtracting the results gives 

0 = </> cos 0 sin 0 sin 0(/j - 1 2 ) + 0(7^ cos 2 0 + / 2 sin 2 0) . (A- 14) 

If the spacecraft is assumed to be dynamically balanced about the z-axis, then I x = / 2 . Making 
this assumption and substituting in Equations A-l 3 and A-l 4, we obtain 

7=^0 (A-l 5) 

and 

I x 6 = 0 . (A-l 6) 

Since I x ¥* 0, 0 must be zero; it follows that 0 is independent of time. Substitution of Equations A-l 5 
and A- 16 in Equation A- 12 yields 


p z = I x <t> cos 0 = / 3 (0 cos 0 + 0) . 

Solving for 0, the angular velocity of the spacecraft about the spacecraft z-axis, we obtain 


(A- 17) 


(/, - /,)0 cos 0 

0= ~7 — • (A- 18) 

2 3 

Solving for 0, the angular velocity of the spacecraft z-axis about the momentum vector or, more 
precisely, the angular velocity of the line of nodes, we obtain 


: = 

* - -7 3 ) cos 0 ' 


(A- 19) 


Let us refer to 0 as precession rate and 0 as the spin rate (to be distinguished from the apparent 
rotation rate of the spacecraft with respect to a fixed external point). (Note that by Equation A-9, 
cu z = 0 when 0 = 0.) This apparent rotation rate has an average value of 0 + 0. The precession cone 
half angle is 0, and the rate of change of the precession cone half angle is 0. If /j// 3 < 1, 9 tends 
toward zero, and the precession coning will damp out in time. Since it is usually desirable for the sat- 
ellite to rotate about the z-axis, most spin-stabilized satellites are balanced so that the z-axis coincides 
with the largest moment of inertia. The rest of this text will deal only with this case, resulting in the 
fact that 0 = 0. It should not be assumed that zero precession implies that 0 goes to zero. Equation 
A-l 9 states that when 0 = 0, then 0 = I 3 i>/(I X ~ / 3 ). 


Equations Determining Spacecraft Aspect 

Figure A-3 shows the spacecraft momentum vector relative to the sun and zenith vector on the 
celestial sphere (which is centered at the spacecraft’s center of mass). The zenith vector is defined as a 
vector from the earth’s center through the spacecraft’s center of mass. The earth’s center is located on 
the celestial sphere at (RA Z + 180) and (~D Z ), where RA z is the right ascension and D z is the declination 
of the zenith vector. The great-circle arc from the sun to the zenith vector is r?. Applying the laws of 
sines and cosines to the spherical triangles in Figure 3, we write the equations for the values of r? and 
the angles L 1 and ZL3 as 
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CELESTIAL NORTH 



Figure A-3— Spacecraft momentum vector relative to the sun and the zenith 

vector. 

cos 77 = sin sin D c + cos D _ cos Z> cos (RA V - RA ) , if 0 < 77 < 1 80° , 

z s z s s z 

cos j3 - cos n cos 5 

cos L 1 = : 7 — 7 , 

sin 77 sin 0 

sinj3sinZ.2 


sin D s - cos 77 sin D z 
sin 7? cos D z ’ 

sin (jRv4 s - RA Z ) cos D s 
sin 77 

where j3 is the angle between the sun and the momentum vector (spin axis) and 5 is the angle between 
the zenith and momentum vectors. Using the above values of L 1 and LI, we can derive the right as- 
cension RA l and declination D L of the momentum vector: 


cos L3 = 

and 

sin L3 = 
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and 


sin D l = sin D z cos 5 + cos D z sin 5 cos (Z.3 + L 1) , if 90° < D L < 90° , 


cos (RA l ~ RA z ) = 
sin ( RA l - RA Z ) = 


cos 5 - sin D z sin D L 
cos D z cos D l 

sin (L 1 + A3) sin 5 
cos D l ’ 


RA l = RA Z + (RA l - RA Z ) . 

Therefore, to derive the right ascension and declination of the spacecraft’s momentum vector, we 
must know several quantities. The right ascensions and declinations of the zenith vector (RA Z ,D Z ) and 
the sun ( RA S ,D S ) are accurately known from earth-based information, whereas )3, 5, and L2 must be 
determined on board the spacecraft. From Figure A-3 it is clear (since apparent rotation rate is 0 + i|/) 
that 

L2 = (0 + 0)A t c ±180°, (A-20) 

where A t c is the time required for the reference plane to move from the sun to the earth’s center. 
Therefore, it is the measurement of (5, 5, and (0 + 0)A t c that the aspect system must perform. 
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Appendix B 


MATHEMATICAL REPRESENTATION OF EULERIAN ROTATIONS 


Any vector N in the X, Y, Z coordinate system, written in column form, may be transformed into 
the £', r)', f' system by the application of the rotation matrix D: 

N x cos 0 sin 0 0 N x 

N n , = D Ny = -sin 0 cos 0 0 N y 

N f , N z 0 0 1 N z 

This vector, when acted upon by the rotation matrix C, yields 

N j N x 1 0 0 cos 0 cos 0 0 N x 

= CD N Y = 0 cos 8 sin 8 -sin 0 cos 0 0 N Y 

N z 0 -sin 8 cos 8 0 0 1 N z 

Finally, application of the rotation matrix B produces the N vector transformed to the x,y, z coordi- 
nate system: 

N x N x cos 0 sin 0 0 1 0 0 cos 0 sin 0 0 N x 

N y = BCD Ny = -sin 0 cos 0 0 0 cos 8 sin 8 -sin 0 cos 0 0 N y 

N z N z 0 0 10 -sin 6 cos 8 0 0 1 N z 

Hence, any vector in the X, Y, Z coordinate system can be determined in the x,y,z coordinate system 

by the application of the transformation matrix A = BCD. Performing the matrix multiplication, we 
obtain 

cos 0 cos 0 - cos 6 sin 0 sin 0 cos 0 sin 0 + cos 8 cos 0 sin 0 sin 0 sin 8 

A = -sin 0 cos 0 - cos 8 sin 0 cos 0 -sin 0 sin 0 + cos 8 cos 0 cos 0 cos 0 sin 8 . (B-l) 

sin 8 sin 0 -sin 8 cos 0 cos 8 
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Since A is an orthogonal transformation, then A~ 1 = A T , and any vector in the x,y,z system can be 
transformed to the X, Y, Z system by application of the transformation matrix A' 1 , we have 


A~ l = 


cos <t> cos \}j - sin <j> cos 0 sin ^ 
sin <p cos + cos 0 cos (j> sin \}/ 
sin 0 sin ^ 


-sin \p cos <f> - cos 0 cos \p sin <j) 
-sin \jj sin 4> + cos 0 cos \}j cos <j> 
sin 0 cos \p 


sin 0 sin <j> 
-sin 0 cos <p . 
cos 0 
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Appendix C 


MOMENTUM AND ANGULAR VELOCITIES 


Momentum 

Consider the momentum vector L of a rotating spacecraft oriented along the Z-axis of the X, Y, Z 
coordinate system. (See Figure A-2a.) Writing the momentum vector as a column vector in the X, Y, 
Z system, we have 



L x 


0 

L = 

l y 

= 

0 


Lz 


L 


Hence, transforming the momentum vector to the x,y, z coordinate system, we obtain the instantane- 
ous values of momentum about the x, y, and z axes. Denoting these values as p x , p y , and p z , respec- 
tively, we have (because A is defined as in Equation B-l) 


Therefore, 


and 


Px 


0 


L sin i// sin 0 

Py 

= AL = A 

0 

= 

L cos sin 0 

Pz 


L 


L cos 0 


p x = L sin \p sin 0 , 
P y = L cos i// sin 0 , 
p z = L cos 0 . 


Angular Velocities 

Consider = d4>ldt = <p as vector directed along the Z-axis (see Figure A-2b). Writing <p as a 
column vector, we obtain 
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Thus, 


and 





d COS 0 


6 cos 0 


- 

II 

•co 

II 

8 sin 0 

= 

-6 sin 0 

co n 


8, 

0 


0 

e z 


z 





co„ = 8 cos 0 , 

a x 

oj„ = -6 sin i p , 
y 

“», = 0 - 
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Finally, consider co^ = d\p/dt = \J/ a vector perpendicular to the plane containing the line of nodes 
and the x-axis and directed along the z-axis. Hence, 




0 


~ 

0 



* 


Now, combining the x, y, and z components of 4>, \p, and 0, we establish a relationship between the 
angular velocities oj x , co y , and co z and the rate of change of the Euler angles: 

co = 0 sin i// sin 0 + 0 cos ^ , 


and 


co y = <j> cos \p sin 0 - 0 sin \j / , 
co = 0 cos 6+ \jj . 
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Appendix D 


CALCULATION OF 5 


The information contained in the relative position of the detector output pulses to the command- 
slit crossing of the sun provides a measure of 5 , the inclination of the spacecraft zenith vector to the 
spin axis. Consider the celestial sphere shown in Figure D-l with the spacecraft located at the center 
of the sphere. The arc length Xj is the smaller of the two possible great-circle arc distances from the 
sun to the horizon. From spherical trigonometry it follows that 

cos X x = cos 0 cos 7 + sin j3 sin 7 cos A . 

There are two values of Xj that satisfy this equation. This ambiguity can be removed if the magnitude 
of angle A is taken into account: 


If n — A > 0 , then 0 < X x < it . 

If 7 r-A< 0 , then 7 r < X 1 < 27r . 

If ir - A = 0 , then X l = 0 + 7 . 

Angle A corresponds to the relative position of the earth’s horizon with respect to the command-slit 
crossing of the solar disk. Hence, X x is uniquely determined, and 

sin X x = y/l~ cos 2 Xj . 

Now, applying the cosine and sine laws of spherical trigonometry to the triangle containing the angle 
X 2 , we have 


and 


cos X 2 


cos 0 - cos X x cos 7 
sin Xj sin 7 


sin X 2 = sin 0 


sin A 
sin Xj ’ 


Therefore, X 2 is also uniquely determined. 


Let r g equal the mean radius of the earth and h be the height of the satellite above the earth (see 
Figure D- 2 c). Then, 

sin p = — ^—r (0 < p < 90°) , 

r„+h 
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MOMENTUM VECTOR 



Figure D-1 —Celestial sphere representation of the single horizon solution for 6 . 


and, from Figure D- 1 , 


cos X 3 


-cos t? - cos X x cos p 
sin X 1 sin p 


Again there are two values of X 3 that satisfy this equation. Also, from Figure D-1, 

X 4 = 2tt - (X 2 ± X 3 ) , 

and X 4 has two possible values corresponding to the two values of X y Also, 

cos 5 = -cos p cos 7 - sin p sin 7 cos X 4 (0 < 5 < 1 80°) , (D-1 ) 

and 5 is doubled-valued. 

The information contained in the spacing of the horizon-detector output pulses and the relative 
position of the earth’s center from the sun provide alternate methods for 5 calculation. 
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Figure D-2— Definition of variables for the two horizon solution for 5. 


Referring to Figure D-2b and considering the plane containing EA and EC, we have 


irAt 


cos- 


a EB 
EA ’ 


(D-2) 


where A t g is the time between horizon pulses, and T E is the rotation period. From Figure D- 2 , various 
other relationships can be derived: 


EB 

- ED 

+ DB , 




(D-3) 

ED 

= OE 

tan (90 - 5) , 

where 

OE = 

OA cos 7 = OH cos y > 

(D-4) 

DB 

cos 

OG 

(90-5) ’ 

where 

OG = 

OH cos p , 

(D-5) 

EA 

= OA 

sin 7 = OH sin 

7 > 



(D-6) 

ED 

= OH 

cos 7 tan (90 

-s). 
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From Equation D-5, 


DB = 


OH cos p 
cos (90 - 6 ) ' 


Hence, from Equation D-2, 


OH cos p 

EB = OH cos 7 tan (90 - 5) + — — p- 

cos (90 - 5) 


( cos p\ 
cos 7 cot 5 + — — - ) 
sin 5 / 


Thus, from Equation D-2, 


Then, from Equation D- 6 , 


OH[ cos 7 cot 


cos p\ 

S+-r-£) 

sin 8 ) 


IT At 

= EA cos — — 

l E 


Therefore, 


that is, 


Let 


COS P 'K&t a 

OH cos 7 cot 5 + — — - = OH sin 7 cos — — 
sm 0 Tr, 


sin 8 cos 7 cot 8 + cos p = sin 5 sin 7 cos - 


7r Ar 


j, > 

1 E 


cos 7 cos 5 + cos p = sin 5 sin 7 cos - 


7tA/ 


Te ' 


d = cos p , 


and 


Then, 


e = cos 7 , 


7rAt fl 

/ = sin 7 cos -^ 7 — 

*E 


d + e cos 5 = / sin 6 = /-y/l - cos 2 S . 


Squaring both sides, we have 

d 2 + 2 ed cos 6 + e 2 cos 2 8 + f 2 cos 2 5 - f 2 = 0 , 
or 

(e 2 + f 2 ) cos 2 8 + 2 ed cos 5 + (d 2 - f 2 ) = 0 . 


(D-7) 


(D- 8 ) 


(D-9) 


(D-10) 


(D- 11) 


(D-12) 
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Thus, 


cos 5 = 


-ed ± \/c? 2 e 2 - (e 2 + / 2 )(tf 2 - / 2 ) 

e 2 + / 2 


that is. 


cos 5 = 


~de ± fy/e 2 + f 2 ~ d 2 

e 2 T7 


Figure D-2 represents a spacecraft rotating with no precession. An earth-horizon sensor is 
mounted 7 degrees off the spacecraft z-axis so as to sweep over a cone which cuts the earth horizons at 
A and C. The inclination 5 of the momentum vector to the subsatellite zenith vector may be ex- 
pressed in terms of spacecraft elevation 7 and horizon-pulse spacing. To summarize, we have just 
proved that 


where 

d = cos p, 


cos 5 = 


~de ± fy/e 2 + f 2 - d 2 
e 2 +J 2 


if 0 < 5 < 1 80° 


(D-13) 


e = cos 7, 
and 

f = sin 7 cos (irAtjT E ). 

Also, from Figures D-l and D-2, 

p = shf 1 [ r E l{r E + h )], 

A t = horizon pulse spacing, 

and 

T e - rotation period. 

Information giving the value of 5 is also contained in the relative position of the horizon pulses 
to the command-slit crossing of the sun. On the spacecraft, the earth horizon sensor has a pencil field 
of view located 7 degrees from the spacecraft z-axis in the plane of the command slit. If the spacecraft 
is spinning so that 0 = 0, horizon pulses will be symmetrical about the instant when this plane crosses 
the center of the earth (see Figure D-3). The center of the earth is located on the celestial sphere at 
( RA Z + 1 80°) and ( -D ) (see Figure 1 2), where RA z is the right ascension and D z is the declination of 
the subsatellite point zenith vector. The great-circle arc from the sun to the earth’s center on the celes- 
tial sphere is 1 80° - 77. The time midway between horizon pulses minus the time the command slit 
crosses the sun is A t . The relationship between 5 and A t can now be written from the spherical tri- 
angles in Figure D-3. 

The angle between 180° - 5 and j3 is L7, i.e., (0 + i^)At c , as is used in Equation D-14, below. 

Since 0 + \p = o> z = 2ir/T, this angle can be expressed as 2ti A t JT\ this is how it appears in Figure D-3. 
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Figure D-3— Relationship between At c and 6. 


Thus, we have 

cos (180° - 7?) = cos j3 cos (180° - 5) + sin j(3 sin (180° - 5) cos [(0 + 0)A t c ] ; 


therefore, 

cos t? = cos (3 cos 5 - sin jS sin 6 cos [(0 + 0)A/ c ] . 

Solving for 5 gives 

rs ±!Xy/? 

cos 5 = ■, , (0 < 6 < 7r) , 

5 2 +M 2 

where 

r = -cos i?, 

s = -cos j3, 
and 

M = sin j3 cos [(0 + 0)At c ]. 

Also, from Figure A-3, 

cos 17 = sin D z sin D s + cos D z cos D s cos (RA s - RA Z ), if 0 < 77 < 1 80° , 

where 

D = declination of subsatellite zenith, 

Z 

D s = declination of sun, 

RA = right ascension of subsatellite zenith, 
and Z 

RA S = right ascension of sun. 


(D-14) 


(D-15) 


(D-16) 


(D-17) 
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Appendix E 


SCHEMATIC DIAGRAMS FOR ASPECT SYSTEM 
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Figure E-3-IMP optical aspect: timing 16-183. 
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Appendix F 


SCHEMATIC DIAGRAMS FOR SSC SYSTEM 
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Figure F-3— IMP SSC, second section. 


















Figure F-6-1MP optical aspect: SSC 16-150. 

























































